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Abstract 

The existence and properties of coherent pattern in the multisoliton solutions of the dKP 
equation over a finite field is investigated. To that end, starting with an algebro-geometric 
construction over a finite field, we derive a "travelling wave" formula for iV-soliton solutions 
in a finite field. However, despite it having a form similar to its analogue in the complex field 
case, the finite field solutions produce patterns essentially different from those of classical 
. interacting solitons 
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Ch ■ 1 Introduction 

> 

There are many diverse methods in the modern theory of integrable systems. One of them, an 
algebro-geometric approach [HJ |2j , was applied to obtain solutions of discrete soliton equations 
over finite fields 01 Q • Within this approach efficient tools for finding algebro-geometric 
solutions based on hyper elliptic curves of arbitrary genus were proposed [HE]. These results are 
in direct analogy to the complex field case, but there are some peculiarities. For instance, since 
finite fields are never algebraically closed there are many more possibilities for the construction 
of breather type solutions (for short discussion see [2]). Also other properties, such as fmiteness 
or cyclic structure, reflect in the character of solutions. 

Our aim here is to discuss the appearance of stable travelling patterns for a general iV-sofiton 
solution over a finite field. To do that we write a general determinant formula (JSJ for iV-soliton 
solution in a travelling wave form (|12|) . The formula obtained is analogous to the typical form of 
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a soliton solution of a Hirota bilinear equation (see e.g. page 23). Since the vacuum solution 
in this setting is r = 1, it is necessary to make a slight modification of the algebro- geometric 
construction presented in j^j. 

The determinant form of the solutions of the dKP equation was already investigated |13| . 
Since the determinant formula is neither a Casorati nor a discrete Gram type determinant, 
we provide a direct link to the travelling wave form without referring to the previous work. 
Moreover, even though we are mainly concerned with finite fields, we point out the calculations 
performed are valid for arbitrary fields. 

There is one more reason to transform the finite field solutions into travelling wave form. 
There exists a well established systematic procedure for deriving soliton cellular automata start- 
ing from discrete soliton equations in Hirota form [141 HUj . Since solutions over a finite field could 
be interpreted also as cellular automata, we need a convenient form to investigate relationships 
between them. 

This article is organized as follows. In Section [21 we recall and transpose some results from 
an algebro-geometric construction of finite field valued solutions of the discrete KP equation 
with explicit determinant formula for iV-soliton solutions 3 to the case of dKP with arbitrary 
coefficients (|2"|). In the next section we prove that the iV-soliton solution can be rewritten 
in travelling wave form. As a remark we also give an alternative proof which explains why 
only pairwise interaction terms appear in the iV-soliton solutions which appear as substitution 
operations for matrix elements. In Section 0] we discuss the patterns produced by finite fields 
solitons and give some examples. The main conclusion is that travelling wave patterns in N- 
soliton solutions obtained by the algebro-geometric approach are generally absent for N > 2. 

2 Algebro-geometric construction of solutions of dKP equation 
over finite fields with simple vacuum. 

A finite field version of an algebro-geometric construction of solutions r : Z 3 — ► F for the discrete 
KP equation 

(Tit)(T 2 3t) - (T 2 t)(T 13 t) + (T 3 t)(T 12 t) = (1) 

was studied in By T, we denoted here a shift operator in a variable rii, for example 
T 2 t(tii, n 2 , n 3 ) = T(ni,n 2 + 1,713). If we prefer to have vacuum solution f(ni, 712,713) = 1 
then we need to consider the dKP equation with coefficients Z% G F, i.e. 

^i(7ir)(T 23 f) - Z 2 {T 2 f)(T l3 T) + Z 3 (T 3 f){T l2 f) = (2) 

with the constraint 

Z 1 -Z 2 + Z 3 = 0. (3) 
Note that a correspondence between r and f is achieved by taking 

r = (Z 3 /d)- nin2 (Z 1 /d)- n2n3 (Z 2 /5)- nin3 T (4) 

for nonzero constants Z±, Z 2 , Z 3 and any nonzero 5 £ F. 

It follows from ;3,, in the case of purely soliton solutions (i.e. from a curve of genus g = 0), 
that A^-soliton solutions can be expressed in a determinant form (Theorem Q below). As a 
component of this formula we need a vacuum solution of Q 

r = (A 1 - A 2 ) n ™{A 1 - A 3 ) n ^(A 2 - A 3 ) n2n \ (5) 
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where Ai £ F, and auxiliary functions </>°, a = 1, 2, ...iV, in the form 



i — C a \ C a — Ai 

k=l 

Parameters C a , where a = 1, . . . ,N, may take value in some finite algebraic extension LdK, 
but they are constrained to the IK-rationality conditions Vcr G £r(L/K), o~(C a ) = C a i. By 
G(L/K) we denoted the Galois group, i.e. the group of automorphisms of L, with fixed field 
K. Similarly, N pairs Dp,Ep € L, for (3 = 1, . . . , N, satisfy the K-rationality conditions Vcr € 
G(L/K.) : a({Dp, Ep}) = {D^ , Ep/}. These conditions give rise to some generalisation of 
breather type solutions (see [£]). We assume all parameters in the construction are distinct. 
Finally, denote by 4>a(D, E) the N x N matrix with element in row [5 and column a given by 

[<j> A (D, E)) aP = 4> a (Dp) - ME/3). (7) 

Theorem 1. The function r(ni, 7t2, 743) given by 

T = T - 1 det^ A (D,E), (8) 

where 7 G F is some constant, is the F-valued N-soliton solution of the discrete KP equation 
©• 

In the case of dKP equation P|). it follows from Q) that a vacuum function fo can be chosen 
in the form 

Mi - A 2 \ m Mi - A 3 \""" / As - .4s V 2 "' 

Then for parameters A\ — A2 = Z3, A\ — A^ = Z 2 , A2 — A3 = Z\ we have f = 1 and 

f = det ^ A {D,E), (9) 
is the F-valued iV-soliton solution of the equation (J2J). 

Remark. The same result can also be derived using the general algebro-geometric construction 
of solutions of the version of the dKP equation given in @. To do this, the wave function ip is 
given by the Definition 1 in [3], but with a different definition of the expansions of ^(ni, ri2, 713) 
at Ai, namely, for i,j, k = 1, 2, 3 with i ^ j ^ k ^ i, 



ip(ni, n 2 , n 3 ) = C J2 ZfZ^lfin^, m)t\ 



s=0 

where U are fixed IK-rational local parameters at Ai. The linear equation for tp in the general 
case is of the form _ 

(Tii;-T j ij) + Z k ^ r ^ = 0. 

Co 

The remaining part of the construction is the same as [3] ■ 
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3 Travelling waves form for the 7V-soliton solution 

In this section we wish to transform the soliton solutions given in (JSJ) into an equivalent but 
more convenient form. In doing this we will use the fact, referred as a gauge invariance, that 
for any solution r of the dKP equation, r' = a ni f3 n2 j n3 S ■ r is also a solution for any constant 
a, /?, 7, 5. We write r' ~ r if r' can be obtained from r using this gauge invariance. 

Lemma 1. Let M(x, y), where x = (x\, . . . , x n ) and y = (yi, . . . , y n ) denote the Cauchy matrix, 
the matrix with (i,j)th entry — yj). It is well known that 



(x p - x q )(y q - y p ) 



Up,a( X P ~ Vq) 



n 



q ( x p - y q ){x q - y P ) Y (x p - y p ) 



(10) 



Also we have, 



M(x, y) _1 M(x, z) = diag 



U p ( x p 



Up^iiVp-Vi) 



M(y,z) diag 



Up(yp-Zi) \ 

U p ( x p -Zi)j ' 



(11) 



Proof. Let xr = (x\, . . . , Xk+i ■ ■ ■ , x n ). Then the (i, j)th entry in the inverse of M(x, y) is 



[M(x,y)" 



U.J 



i+? . detM(x 3 -,y t - 
1 J det M(x, y) 



p+i 



( X t 



n 



y P ) 



Further, the (i,j)th entry in the product M(x, y) 1 M(x, z) is 



Up^iivp -yi)f^( x k- zj) U P M Xk ~ 
U p ( x p - yi) 



^(-ir- fc y(x^n(^-yp)II( 3 

p-^-i P~£k 



where V(x) = Y\ p>g (x p — x q ) is the Vandermonde determinant in variables xx,...,x n and V^(x^) 
is the Vandermonde determinant in variables x\, . . . , Xk-i,Xk+i, ■ ■ ■ , x n . In the (i, j)th entry, the 
sum is of degree (n — l)(n — 2)/2 + 2(n — 1) = n{n — l)/2 + n — 1 and we will factorize it by 
identifying all of its zeros. First, there are n(n — l)/2 coming when x r = x s for any r < s and 
n — 1 from y r = Zj for r ^ i. This argument uses the fact that for any polynomial function / 
of degree less than n — 1, one has Ylk=i(~ l) n_fc ^'( x fc)/( a; fc) = since the left hand side is the 
expansion by the final column of the vanishing determinant 



1 X\ 



1 x r 



„n-2 



„n-2 



f{ x r 



Hence the (i,j')th entry is, up to a constant factor, 

U P ( x P-yi) Up^i(yp-zj) _ U p ( x p 



yi) i U P (y P - zj) 



Up^iivp-yi) U P ( x P-zj) Up^i(y P - Vi) Vi - ZjU P ( x P- zj)' 

To verify that it is precisely correct, and hence complete the proof of (|llj) . we observe that, as 
it should, the above equals Sy, the Kronecker delta, when z = y. □ 
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Theorem 2. Let q denote any fixed generator o/F* i.e. a multiplicative subgroup of the finite 
fieldW. The N-soliton solution © of the dKP equation over a finite field ¥ admits the following 
form 



r'= £ (-1) #J II a* (12) 

JC{1,...,N} \i,i'eJ; i<i' J 

where the sum is taken over all subsets of {1, . . . , N} and j^J denotes the cardinality of J. In 

(A - Dj)(Ei - Ej) 



lJ ' (A - A) (A - EiY 



(13) 



J 3 ) 

the exponents are fjj = rjj + rjj where 



j 

3 



fc=l 

and the parameters p k ; and phase constants rff are defined by 



N j—. •, AT 



q -Di-A* md Q - 11 ( Cp -Ek) 11 ,(VAj' (15) 

The name travelling wave form comes from the form of the q^ , which are in direct analogy 
with the usual term exp(kx — tot) of linear plane waves. Note, that for any r satisfying (J2J), 
t(1, 0, 0), r(0, 1,0) and r(0,0, 1) depend only on cross ratio of appropriate points on the pro- 
jective line. We also point out that to find p\ and 77° ■ we need to solve a discrete logarithm 
problem. 

Proof. The determinant in Q is 

r = (-m(c, d , ^ (n (|^)"') + m(c, e) ^ (n (t^f)) ■ 

where M denotes the Cauchy matrix as defined in Lemma ^ and so 



f ~ T > ■- det (i - M(C, D)- l M(C, E) diag m ^ 



A-^ N "* 



.A - -4a 

Using and the fact that for any matrices P, Q, det PQ = det QP, we get 

t' = det (J - M(D,E) diag(A - A) diag(^))- 

Now the determinant det(P + Q), of the sum of two N x N matrices, may be expressed as 
the sum over all subsets J C {1, . . . ,N} of detPj where the ith column of Rj equals the zth 
column of Q if i £ J and otherwise equals the ith column of P. Using this fact, 

T '= E (-l^detM^A^II^-WS 

JC{1,-,N} ieJ 

where Mj denotes the Cauchy matrix with row and column indices restricted to J. The formula 
(j!2|) now follows immediately from (jlOj) . □ 
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Remark. We also present here an alternative proof of Theorem |21 The argument presented 
shows why there is only pairwise interaction of solitons. 

Starting from for each i one divides the ith column of by (pj(Di) and we see that 



t ~ t = det where 



i.i 



Di-Cj 

Di - a 



Di-Cj 
Ea — d 



n 

k=l 



Dj - A k 
Di - Ai 



in: 



where r]j is given by (|14() . is also solution of ©. Further, if we define 

D{ — Ci\ I Di — Ci 



ij 



D< - d 



Ej — d 



(16) 



(17) 



it is easy to see that det <£> ~ det <fe. 

We will next show how det $ can be obtained from a vacuum solution det where [<&o]?j = 

( d'-c' ) means of substitution operations. The key idea we use in finding coefficients of 
different powers of q in the formula (|12|) is observation that, according to (|17|) . for each term 
including q Vk we should make the substitution 



Di-Cj 
Dh — Ci 



Di-Cj 
Eh — Ci 



(18) 



for each i, in the kth row of $o- Then det $ can be expressed in terms of determinants obtained 
by making appropriate replacements in det<l>o- 
First notice that 

det <f>o = M(-C, -D) JJ(A - Cj) = f[ ^zM^lM . (19) 



i,i=i; i<j 



Denote by $j the matrix with replacements (|18jl in the rows A; € J. We first consider the 
case J = {k} for fixed k. Since Dk appears only in the kth row and kth column, then $j after 



multiplying kth column by 
Then 



-c k 



D k -C k 



becomes the determinant (fT§|) with Dk replaced by E^. 



det$ 



Dk — Cfc 



Ek 


— Ck 


D k 


— Ck 


Ek 


— Ck 


D k 


- Ck 


Ek 


— Ck 




det$ 



ci- 



In the notation of the theorem, the coefficient of is det &{k} = —q Vk det $o- 

Now we consider coefficient of q^k^jVk _ This is the determinant det $ j obtained from det $0 
by making replacements (|18|) in rows k € J for general J C I. The number of elements in J is 
denoted by # J. Repeating the procedure described above for each k £ J we arrive at 



det$ 



n 

vfceJ 



Ek — Ck 



det$ 



D k ^E k ; k£j 



(20) 
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Then we have 



det$ 



(D kl - Dk 2 )(E kl - E k2 

D k ^E k ; fcgj \ ^-L (E kl — Dk 2 )(Dk 1 — E k2 

\Ki ,Kj f(ZJ 



n 




det$ 



o- 



The extra terms for pairs {ki,kj} £ J are present because, since kj ^ ki, we also need to 

perform the replacements D k . — ► Ej^ in the factor n^=i- j'^kt ( 5^ -c'', ) ( -a' ) • This 

is done by multiplying by \^ kl _n k2 \m kl -^° 2 \ ■ We point out here the extra corrections are for 
pairs {/cj, kj} £ J and it is this that is responsible for the existence of only pairwise interaction 
terms in the formula (|12|) . So finally, the coefficient of q^jeJ^i is 

det^ = (-1)* J I J] | [HqA det4» . (21) 

\i,i'eJ; »<i' / VfceJ / 

After dividing all terms by det $o an d collecting them together we obtain formula (1121) . 

Notice that our considerations are valid for any field of definition of the parameters C, and 
Dj. 

Remark. The results obtained above apply also to the Discrete Analogue of Generalized Toda 
Equation jS] 

^ , lfmi+l,m2,mzfmi — l,m2,rnz ^2fm-i,m2+l,m^fmi,rn2 — l,mz ^3/mi,m2,m3+l/mi,m2,m3-l — 0. 

The transition from dKP to DAGTE (where we assume (mi + m<i + 711,3) = mod 2) is given 
by f mi ,m2,m3 ■= r(m, n 2j n 3 ) where = \{rrii - mj - m k ) for i / j / k / i, or equivalently 
rrii = —(rij + n k ) for i ^ j ^ k ^ i. In this situation /o,o,o = t(0, 0, 0). With the variables 
1111,7112,1713 the exponent of a travelling wave given by Q14[) have similar form, namely 77^ = 
Y^j=iPi rn j > wn ere = ^fc=i(^^jfc — and is the Kronecker delta. 



4 Pattern structures in soliton interactions 

So far we have seen that the finite field iV-soliton solutions we obtained have the same structure as 
the complex field case. In particular, the travelling wave form (|12[) is completely analogous to the 
complex field case. Despite this correspondence, the typical soliton-like interaction pattern is not 
present in the finite field case. Below we argue that one cannot expect that an arbitrary iV-soliton 
solution over a finite field will be a collection of asymptotically separated waves interacting in 
the way characteristic of the complex field case. 

These differences in the interaction properties follow naturally from the differences in struc- 
ture between finite fields and the complex numbers. First, finite fields have no good (total) 
ordering which are consistent for both addition and multiplication and so it is impossible to find 
in the finite filed case a direct analogue of wave amplitude; rather than follow a wave propagation 
by observing how its points of maximum amplitude move we may trace only the propagation of 



7 



patterns. The appearance of a one-soliton solution however in the finite field case mimics the 
usual appearance quite well, as can be seen on Figure ^ 

The next difference results from the replacement of exponentials in the complex case by 
powers of a generator q of the multiplicative group F* of the finite field F. Since g' F ' _1 = 1, this 
implies periodicity of r(ni, 712,^3) with respect to each variable rtj. As a consequence we can 
not discuss the asymptotic behaviour of the solution in usual sense. Moreover, for any given 
solution one could restrict analysis to a finite base cube containing all information about the 
solution, the rest is but periodic repetition. The length of any edge of a base cube is at most 
|F| - 1. 

Considering Theorem|21 the ith one-soliton component of the iV-soliton solution is unchanged 
by a shift in the lattice by n l = (n\, n 2 , n 3 ) for any n\, n 2 , n 3 satisfying = 1, or equivalently, 

3 

Vi = J2Pi n k = ° mod(|F|-l). (22) 

k=l 

Since expression (|12j) contains q m for i € {1, 2, . . . , N}, a period vector n = (rii, ?i2, 713) for this 
solution should be a common solution of (|22|) for all i. In general, it is impossible to find a 
nonzero solution for N > 3 and it means there is no additional structure within the base cube 
in this case. An example of such a 3-soliton solution is shown in Figure |3 and the structure of 
a base cube is discussed below in more detail. 

Examples. Fix the finite field to be F = F17, i.e. the field of integers modulo 17. As a 
generator of F* we choose q = 3. Let us fix A± = 7, A2 = 4 and ^3 = 3 so that the coefficients 
in the dKP equation are A2 — A3 = Z\ = 1, A\ — A 3 = Z<i = 4 and A± — A2 = Z 3 = 3. 

One-soliton solutions. In Figure^we present three one-soliton solutions of the dKP equation 
over F. For the solution A on the left, denoted by subscript 1 we have chosen C\ = 11, D\ = 6, 
E\ = 9; for the solution B in the middle, denoted by 2, we have C2 = 10, D2 = 12, E2 = 14 and 
finally C3 = 8, D 3 = 13, E3 = 15 for the solution C on the right. The respective parameters in 
the exponent ()14|) are 

Pi = (p\,pIpD = (6,7,14), p 2 = (pipl P 3 2 ) = (6,9,5), p 3 = {p\,pl,pl) = (11,5,10). 

Thus periods in variables ni,U2,ns are: 8,16,8 for the soliton 1, 8,16,16 for the soliton 2 and 
16, 16, 8 for 3. We fixed n\ = 0, since increasing any variable by 1 results in shift in the other two 
variables and so the plots for any value of n\ is simply a translation of the plots we present. This 
comes from the fact that in general there exists a nonzero solution in two variables nj,nk of the 
single equation (|22|) with the third variable rej being fixed (i, j, k € {1, 2, 3}, i ^ j ^ k ^ i). The 
special case is ifp 1 is not a linear combination of pi andp k with coefficients from Z mod (|F| — 1). 
(This could happen for instance if p> and p k are zero divisors of |F| — 1.) In the case presented 
in Figure E period vectors might be chosen to be n\ = (1,0,3) for A, H2 = (1,0,2) for B and 
n3 = (1,1,0) for C. Further, the period vectors n for these solutions in the plane n2,n% are 
n\ a = n 3a = (0,2,-1), n lb = n 3b = (0,0,8), n 2a = (0,3,1) and n 2 b = (0,-1,5). (All periods 
vectors are not uniquely determined.) 

Two-soliton solutions. In Figure [2] the two-soliton solutions representing the interaction of 
two of the one-soliton solutions shown in Figure ^ are presented. The upper line is for n± = 
while the lower is for n\ = 1. In the cases A interacting with B (called AB) and B interacting 
with C (AC) it is a simple shift by respective vectors (1, —1, —4) and (1, —5, —5). (Note: Since 
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Elements of F17 are represented on the scale below: from - dark to 16 - light gray. 




2.5 5 7.5 10 12.5 15 2.5 5 7.5 10 12.5 15 2.5 5 7.5 10 12.5 15 



Figure 1: A plot of t (711,122,713) function of a three one-soliton solutions (see Examples for 
details). We fix n\ = 0, and 112,113 G {0, . . . , 16}. The 722 axis is directed to the right and the 
723 axis is directed upwards. 




2.5 5 7.5 10 12.5 15 2.5 5 7.5 10 12.5 15 2.5 5 7.5 10 12.5 15 



Figure 2: A plot of r(ni, 712,723) for the three two-soliton interactions (AB, AC, BC) of one- 
soliton solutions presented in Figure Q We fix n\ = 0, and 722,713 G {0, . . . , 16}. In the second 
row we have m = 1. The tt-2 axis is directed to the right and the 723 axis is directed upwards. 
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Elements of F17 are represented on the scale below: from - dark to 16 - light gray. 




2.5 5 7.5 10 12.5 15 2.5 5 7.5 10 12.5 15 



Figure 3: A three-soliton solution r (711,712,^13) being the solitonic sum of those from Figure ^ 
for m = 0, 1, 2,4 and 8. Axes: ri2 directed to the right, 713 directed upward. 

g 16 = 1, the lines n = and n = 16 are identical but both of them are shown on plots.) For the 
case AC there are no such a shift, because the two equations ()22l) for i = 1 and i = 3 for n<i and 
ri3 has no nonzero solution for any value of n\. This is because n\ ■ (6, 11) Zi6 • (7, 5). Since 
span(ni a , nib) H span(n*2a, n2b) = {0} there are no nonzero period vectors in the plane ri2, ri3 for 
either AB or BC. On the contrary, for AC the period vectors in this plane are exactly the same 
as for soliton A and C (since they are the same). Similarly one could examine other planes 
obtaining plane period vectors (0,0,8) for AC, (8,0,0) and (4,8,0) for AB and none for BC. 
Periods for AB are (8,16,16), for AC are (16,16,8) and for BC are the maximal (16,16,16). 
Notice that BC shows that it is possible for a two-soliton solution to have no structure within 
the base cube. 

A three-soliton solution. In this case there are three equations ()22|) for three variables, so in 
general there are no nonzero solutions. The example presented in Figure 01 using parameters in 
A, B and C is of this kind. Because of this, there are no extra period vectors and so this solution 
has no structure within a base cube. 

It is clear that this is also the generic situation for the interaction of A r -soliton solutions 
obtained by the algebro- geometric approach for N > 2. In short, we have seen that the inter- 
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action of three or more finite field soliton solutions of the dKP equation has no more structure 
in the base cube than random data. Given a finite field F of sufficient size, it would seem to be 
computationally impractical to attempt to reconstruct the parameters C a , D a and E a which 
define the solution, from the seemingly random data that they generate. This leads one to 
imagine possible applications of such solutions in encrypted data transmission. In conclusion, 
we note that applications of elliptic curves over finite fields have already been studied for some 
time CU 
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